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Different future long baseline neutrino oscillation setups are discussed and the remarkable potential for very 
precise measurements of mass splittings and mixing angles is shown. Furthermore it will be possible to make 
precise tests of MSW effects, which allow to determine the sign of Am^. Finally strong limits or measurements 
of leptonic CP violation will be possible. 



1. Introduction 

The existing evidence for atmospheric neutrino 
oscillations includes some sensitivity to the char- 
acteristic L/E dependence of oscillations and 
there is no doubt that the observed flavour transi- 
tions are due to neutrino oscillations. Solar neu- 
trinos were also shown to undergo flavour tran- 
sitions H,^. This solves the long standing solar 
neutrino problem, even though the characteris- 
tic L/E dependence of oscillation is in this case 
not yet established. However, oscillation is un- 
der all alternatives by far the most plausible ex- 
planation and global oscillation fits clearly favour 
the so-called LMA solution for the mass splittings 
and mixings The CHOOZ reactor exper- 

iment provides moreover currently the most 
stringent upper bound for the sub-leading Ue3 el- 
ement of the neutrino mixing matrix. The global 
pattern of neutrino oscillation parameters seems 
therefore quite well known and one may ask how 
precise future experiments will ultimately be able 
to measure mass splittings and mixings and what 
can be learned from such precise measurements. 

The characteristic length scale L of oscilla- 
tions is given by Am^L/Ei, = 7r/2. The atmo- 
spheric Amgj^-value leads thus to an oscillation 
length scale Latm as a function of energy. For 
Arn^i ~ 3 • 10^'' eV, and for neutrino energies of 
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E^ ~ 10 GeV, one finds Latm ^ 0(2000) km, 
i.e. distances and energies which can be real- 
ized by sending neutrino beams from one point 
on the Earth to another. The solar Am|]^ is for 
the favoured LMA solution about two orders of 
magnitude smaller than the atmospheric Am^i, 
resulting for the same energies in oscillations at 
scales Lsoi — (10— 1000) - Latm- Solar oscillations 
will thus not fully develop in LBL experiments 
on Earth, but sub-leading effects play neverthe- 
less an important role in precision experiments. 
Another modification comes from the fact that 
the neutrino beams of LBL experiments traverse 
the matter of the Earth. Coherent forward scat- 
tering in matter must therefore to be taken into 
account, which makes the analysis more involved, 
but it offers also unique opportunities. 

The existing K2K experiment as well as MI- 
NOS and CNGS, which are both under construc- 
tion, are a promising first generation of LBL ex- 
periments. We will discuss here (based on [|j) 
the remarkable potential byond this first genera- 
tion. One important point is that the increased 
precision will allow to test in detail the three- 
flavouredness of oscillations. We will also see 
that it will be possible to limit or measure 9i3 
drastically better than today. Next it will be 
possible to study in detail MSW matter effects 
and to extract in this way sign{Am'^i), i.e. the 
mass ordering of the neutrino states. For the now 
favoured LMA solution it will be possible to mea- 
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sure leptonic CP violation [0. The precise neu- 
trino masses, mixings and CP phases which can 
be obtained in this way are extremely valuable 
information about flavour physics, since unlike 
for quarks these parameters are not obscured by 
hadronic uncertainties. These parameters can be 
evolved with the renormalization group e.g. to the 
GUT scale, where the rather precisely known pa- 
rameters can be compared with mass models neu- 
trino masses and mixings. Leptonic CP violation 
is moreover related to leptogenesis, the currently 
most plausible mechanism for the generation of 
the baryon asymmetry of the universe. LBL ex- 
periments offer therefore in a unique way access 
to precise knowledge on extremely interesting and 
valuable physics parameters. 

2. Beams and Detectors 

LBL experiments have the advantage that both 
source and detector can be kept under precise 
conditions. This includes amongst others for the 
source a precise knowledge of the mean neutrino 
energy E^, the neutrino flux and spectrum, as 
well as the flavour composition and contamina- 
tion of the beam. Another important aspect 
is whether neutrino and anti-neutrino data can 
be obtained symmetrically such that systematical 
uncertainties cancel in an analysis. Precise mea- 
surements require also a sufficient luminosity and 
a detector such that enough statistics can be ob- 
tained. On the detector side one must include fur- 
ther issues, like the detection threshold function, 
energy calibration, resolution, particle identifi- 
cation capabilities (flavour, charge, event recon- 
struction, understanding backgrounds). Another 
source of uncertainty in the detection process is 
the knowledge of neutrino cross-sections, espe- 
cially at low energies Source and detector 
combinations of a future LBL experiment are fur- 
thermore constraint by the available technology 
and one should keep potential improvements of 
source and detector developments in mind. An 
example is given by liquid argon detectors like 
ICARUS ||l2| which are not included in this study, 
but might become extremely valuable detectors or 
detector components in the future, improving the 
potential even further. 
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Figure 1. Neutrino production, oscillation and 
detection via charged current interactions for con- 
ventional beams and superbeams. Most interest- 
ing are the i^^ disappearance and — > Vg 
appearance channels. The i^e beam contamina- 
tion at the level of < 1% limits the ability to de- 
termine the — > Ve appearance oscillation, since 
it produces also electrons. The channel is not 
shown here, but it would become very important 
if tau lepton detection were feasible. 



The first type of considered sources are con- 
ventional neutrino and anti-neutrino beams. An 
intense proton beam is typically directed onto a 
massive target producing mostly pions and some 
K mesons, which are captured by an optical sys- 
tem of magnets in order to obtain a beam. The 
pions (and K mesons) decay in a decay pipe, 
yielding essentially a muon neutrino beam which 
can undergo oscillations as shown in fig. ^ Most 
interesting are the i'^ — s- disappearance chan- 
nel and the i'^ Ve appearance channels. The 
neutrino beam is, however, contaminated by ap- 
proximately 0.5% electron neutrinos, which also 
produce electron reactions in the disappearance 
channel, limiting thus the precision in the ex- 
traction of — > Vf, oscillation parameters. The 
energy spectrum of the muon beam can be con- 
trolled over a wide range: it depends on the in- 
cident proton energy, the optical system, and the 
precise direction of the beam axis compared to 
the direction of the detector. It is possible to 
produce broad band high energy beams, such as 
the CNGS beam or narrow band lower 

energy beams, such as in some configurations of 
the NuMI beam |]l5j . Reversing the electrical cur- 
rent in the lens system results in an anti-neutrino 
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beam. The neutrino and anti-neutrino beams 
have significant differences such that errors do not 
cancel systematically in ratios or differences. The 
neutrino and anti-neutrino beams must therefore 
more or less be considered as independent sources 
with different systematical errors. 

So-called "superbeams" are based on the same 
beam dump techniques for producing neutrino 
beams, but at much larger luminosities ]l3| , |l^ , p^ 
|6j . Superbeams are thus a technological extrap- 
olation of conventional beams, but use a proton 
beam intensity close to the mechanical stability 
limit of the target at a typical thermal power of 
0.7 MW to 4MW. The much higher neutrino lu- 
minosity allows the use of the decay kinematics 
of pions to produce so-called "off-axis beams" , 
where the detector is located some degrees off 
the beam axis. This reduces the neutrino flux 
and the average neutrino energy, but leads to 
a more mono-energetic beam and a significant 
suppression of the electron neutrino contamina- 
tion. Several off-axis superbeams with energies 
of about 1 GeV to 2 GeV have been proposed in 
Japan America |§, and Europe ^d^. 

The most sensitive neutrino oscillation chan- 
nel for sub-leading oscillation parameters is the 
Vfj, — > appearance transition. Therefore the de- 
tector should have excellent electron and muon 
charged current identification capabilities. In 
addition, an efficient rejection of neutral cur- 
rent events is required, because the neutral cur- 
rent interaction mode is flavor blind. With low 
statistics, the magnitude of the contamination it- 
self limits the sensitivity to the v^, transi- 
tion severely, while the insufficient knowledge of 
its magnitude constrains the sensitivity for high 
statistics. A near detector allows a substantial 
reduction of the background uncertainties p7| , p2[ 
and plays a crucial role in controlling other sys- 
tematical errors, such as the flux normalization, 
the spectral shape of the beam, and the neutrino 
cross section at low energies. At energies of about 
1 GeV, the dominant charge current interaction 
mode is quasi-elastic scattering, which suggests 
that water Cherenkov detectors are the optimal 
type of detector. At these energies, a baseline 
of about 300 km would be optimal to measure at 
the first maximum of the oscillation. At about 



2 GeV, there is already a considerable contribu- 
tion of inelastic scattering to the charged current 
interactions, which means that it would be useful 
to measure the energy of the hadronic part of the 
cross section. This favors low-Z hadron calorime- 
ters, which also have a factor of ten better neu- 
tral current rejection capability compared to wa- 
ter Cherenkov detectors In this case, the 
optimum baseline is around 600 km. The matter 
effects are expected to be small for these experi- 
ments for two reasons. First of all, an energy of 
about 1 GeV to 2 GeV is small compared to the 
MSW resonance energy of approximately 13 GeV 
in the upper mantle of the Earth. The second 
reason is that the baseline is too short to produce 
significant matter effects. 

The second type of beam considered are so- 
called neutrino factories, where muons are stored 
in the long straight sections of a storage ring. 
The decaying muons produce muon and electron 
anti-neutrinos in equal numbers |^3| . The muons 
are produced by pion decays, where the pions 
are produced by the same technique as for su- 
perbeams. After being collected, they have to be 
cooled and reaccelerated very quickly. This has 
not yet been demonstrated and it is a major tech- 
nological challenge. The spectrum and flavor con- 
tent of the beam are completely characterized by 
the muon decay and are therefore very precisely 
known |2^. The only adjustable parameter is the 
muon energy i?^, which is usually considered in 
the range from 20 to 50 GeV. In a neutrino fac- 
tory it would be possible to produce and store 
anti-muons in order to obtain a CP conjugated 
beam. The symmetric operation of both beams 
leads to the cancellation or significant reduction 
of errors and systematical uncertainties. We will 
discuss in the following the neutrino beam, which 
implies always implicitly - unless otherwise stated 
- the CP conjugate channel. 

The decay of the muons and the relevant os- 
cillation channels are shown in fig. |^. Amongst 
all flavors and interaction types, muon charged 
current events are the easiest to detect. The ap- 
pearance channel with the best sensitivity is thus 
the Ve V transition, which produces so called 
"wrong sign muons" . Therefore, a detector must 
be able to very reliably identify the charge of a 
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Figure 2. Neutrino production, oscillation and 
detection via charged current interactions for a 
neutrino factory for one polarity. 77^ and are 
produced in equal numbers from /i-decays and 
can undergo different oscillations. The —^ 
and channels are most interesting for de- 

tectors with ^ identification. Note, however, that 
excellent charge identification capabilities are re- 
quired to separate "wrong sign muons" and "right 
sign muons". The Vr oscillation channel is not 
shown here, but it would become important for 
detectors with tau identification capabilities. 



muon in order to distinguish wrong sign muons 
in the appearance channel from the higher rate 
of same sign muons in the disappearance chan- 
nels. The dominant charge current interaction 
in the multi-GeV range is deep-inelastic scat- 
tering, making a good energy resolution for the 
hadronic energy deposition necessary. Magne- 
tized iron calorimeters are thus the favored choice 
for neutrino factory detectors. In order to achieve 
the required muon charge separation, it is neces- 
sary to impose a minimum muon energy cut at ap- 
proximately 4 GeV (2^,^ . This leads to a signif- 
icant loss of neutrino events in the range of about 
4 GeV to 20 GeV, which means that a high muon 
energy of Efj_ = 50 GeV is desirable. The first 
oscillation maximum lies then at approximately 
3 000 km. Matter effects are sizable at this base- 
line and energy, and the limited knowledge of the 
Earth's matter density profile becomes an addi- 
tional source of errors. 



3. The Oscillation Pramework 

The general three neutrino oscillation probabil- 
ities in matter are quite lengthy expressions. It is 
useful to simplify them such that an analytic un- 
derstanding becomes possible. Quantitative re- 
sults should, however, be obtained in numerical 
calculations using the full expressions. The key to 
an analytic simplification is to expand the oscilla- 
tion probabilities in small quantities, namely a = 
Am|i/Am§i ~ lO^^ and sin2 26'i.c( < 0.1. The 



matter effects can be parametrized by the dimen- 
sionless quantity A = A/Amli ~ 2VE/A 
where V = \/2GFne- Using A = 
ing terms in this expansion are |2 
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'^ij 1 ^ij — sin '20'ij , 
— i^uo^uij, ss = sins, cs = C0S5, and where in 



where Sij = sin0y, Cy = cos 
Cij = cos 26'; 

eq. (|^) stands for neutrinos and "— " for anti- 
neutrinos. The most important feature of eq. (^) 
is that all interesting effects in the Ve i^^i transi- 
tion depend crucially on 613 . The size of sin^ 2^13 
determines thus if the total transition rate, mat- 
ter effects, effects due to the sign of Am§i and 
CP violating effects are measurable. One of the 
most important questions for future LBL experi- 
ments is therefore how far experiments can push 
the 9i3 limit below the current CHOOZ bound of 
approximately sin^ 2^13 < 0.1. 

Before we discuss further important features of 
eqs. (|l]) and (||) in more detail we would like to 
comment once more on the underlying assump- 
tions and the reliability of these equations. First 
eqs. (|l|) and (^) are an expansion in terms of 
the small quantities a and sin 2^13. Higher order 
terms are suppressed at least by another power 
of one of these small parameters and these cor- 
rections are thus typically at the percent level. 
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The matter corrections in eqs. (|l]) and (H) are 
derived for constant average matter density. Nu- 
merical test have shown that this approximation 
works quite well as long as the matter profile is 
reasonably smooth. A number of very interest- 
ing effects existing in general non-constant matter 
distributions are therefore only small theoretical 
uncertainties. An example is given by asymmet- 
ric matter profiles, which lead to interesting T- 
violating effects |Q , but this does not play a role 
here since the Earth is sufficiently symmetric. 

Note that all results which will be shown later 
are based on numerical simulations of the full 
problem in matter. These results do therefore not 
depend on any approximation. Eqs. (|l|) and (|^) 
will only be used to understand the problem an- 
alytically, which is extremely helpful in order to 
oversee the six (or more) dimensional parameter 
space. The full numerical analysis and eqs. (|^) 
and (H) depend, however, on the assumption of 
a standard three neutrino scenario. It is thus as- 
sumed that the LSND signal |0 will not be con- 
firmed by the MiniBooNE experiment j32). Neu- 
trinos could in principle decay, which would make 
the analysis much more involved. It is assumed 
in this article that neutrinos are stable, and a 
combined treatment of oscillation and decay 
would be much more involved. Neutrinos might 
further have unusual properties and might, for ex- 
ample, violate CPT. In that case neutrinos and 
anti-neutrinos could have different properties and 
LBL experiments can give very interesting limits 
on this possibility , but we will assume in this 
study that CPT is preserved. 

4. Correlations and Degeneracies 

Eqs. (|l|) and (||) exhibit certain parameter cor- 
relations and degeneracies, which play an impor- 
tant role in the analysis of LBL experiments, and 
which would be hard to understand in a purely 
numerical analysis of the high dimensional para- 
meter space. The most important properties are: 

• First we observe that eqs. (|l|) and (||) 

depend only on the product a ■ sin20i2 
or equivalently Am|]^ • sin 2^12. This are 
the parameters related to solar oscillations 
which will be taken as external input. The 



fact that only the product enters, implies 
that it may be better determined than the 
product of the measurements of Am|]^ and 
sin 2^12. 

• Next we observe in eq. (||) that the sec- 
ond and third term contain both a fac- 
tor sin(y4.A), while the last term contains 
a factor sin^(AA). Since AA = 2VL, 
we find that these factors depend only on 
L, resulting in a "magic baseline" when 
2VLmagic = / , whcrc sin(AA) van- 
ishes. At this magic baseline only the first 
term in eq. survives and P{i'e Vfj) 
does no longer depend on (S, a and sin 2^12. 
This is in principle very important, since 
it implies that sin^ 2^13 can be determined 
at the magic baseline from the first term of 
eq. (^) whatever the values and errors of i5, 
a and sin 20i2 are. For the matter density 
of the Earth we find 

L„,ag^c = 7r/4l^ ~ 8100 km , (3) 

which is an amazing number, since the value 
of V could be such that Lmagic is very dif- 
ferent from the scales under discussion. 

• Next we observe that only the second and 
third term of eq. (||) depend on the CP 
phase (5, and both terms contain a factor 
sin 2^13 • a, while the first and fourth term 
of eq. (H) do not depend on the CP phase 
5 and contain factors of sin^ 20i3 and a^, 
respectively. The extraction of CP viola- 
tion is thus always suppressed by the prod- 
uct sin 2013 • ce and the CP violating terms 
are furthermore obscured by large CP in- 
dependent terms if either sin^ 26113 ^ «^ 
or sin^ 2^13 ^ a^. The determination of 
the CP phase S is thus best possible if 
sin^ 26*13 - - a^- 

• The last term in eq. d), which IS propor- 
tional to = (Am2i)^/(Am|i)^, domi- 
nates for tiny sin^ 26*13. The error of Am|i 
limits therefore for small sin^ 26*13 the para- 
meter extraction. 



Eqs. (|l|) and ^ suggest the existence of 
degeneracies, i.e. for given L/E parameter 
sets with identical oscillation probabilities. 
An example is given by a simultaneous re- 
placement of neutrinos by anti-neutrinos 
and Ato§;^ — > — Am^^. This is equivalent 
to changing the sign of the second term of 
eq. (^) and replacing a ~* —a and A 
-A, while A. Eqs. (0) and (||) are un- 
changed, but this does not constitute a de- 
generacy, since neutrinos and anti-neutrinos 
can be distinguished experimentally. 

The first real degeneracy |35j can be seen in 
the disappearance probability eq. , which 
is invariant under 623 — > 7r/2 — 623. Note 
that the second and third term in eq. (|^) are 
not invariant under this transformation, but 
this change in the sub-leading appearance 
probability can approximately be compen- 
sated by small parameter shifts. However, 
the degeneracy can in principle be lifted 
with precision measurements in the disap- 
pearance channels. 

The second degeneracy can be found in the 
appearance probability eq. in the (S — 
0i3)-plane In terms of ^13 (which is 

small) and S the four terms of eq. (||) have 
the structure 

+013 ■ (± sin 6F2 + COS SF3) + F4 (4) 

where the quantities Fi, i = 1,..,4 con- 
tain all the other parameters. The re- 
quirement P{h'e ^ f/^) = const, leads for 
both neutrinos and anti-neutrinos to pa- 
rameter manifolds of degenerate or corre- 
lated solutions. Having both neutrino and 
anti-neutrino beams, the two channels can 
be used independently, which is equivalent 
to considering simultaneously eq. (jj) for 
F2 = and F3 = 0. The requirement 
that these probabilities are now indepen- 
dently constant, i.e. P{i'e i^fi) — const. 
for = and F3 = leads to more con- 
straint manifolds in the ((5 — 0i3)-plane, but 
some degeneracies still survive. 



• The third degeneracy js^ is given by the 
fact that a change in sign of Am|]^ can es- 
sentially be compensated by an offset in 5. 
Therefore we note again that the transfor- 
mation Amgj^ — > — Atojj^ leads to a ^ —a, 
A -A and A -A. All terms of the 
disappearance probability, eq. (|]), are in- 
variant under this transformation. The first 
and fourth term in the appearance proba- 
bility eq. (|l|), which do not depend on the 
CP phase 5, are also invariant. The second 
and third term of eq. (|l|) depend on the 
CP phase and change by the transforma- 
tion Am|]^ — Am|]^. The fact that these 
changes can be compensated by an offset in 
the CP phase 5 is the third degeneracy. 

• Altogether there exists an eight-fold degen- 
eracy as long as only the v^, 
77^ — > T7^ , Vf, ^ I'fi and Fg —^ 17^ channels 
and one fixed L/E are considered. How- 
ever, the structure of eqs. (||) and (||) makes 
clear that the degeneracies can be broken 
by using in a suitable way information from 
different L/E values. This can be achieved 
in total event rates by changing L or E 

but it can in principle also be done 
by using information in the event rate spec- 
trum of a single baseline L, which requires 
detectors with very good energy resolution 
. Another way to break the degeneracies 
is to include further oscillation channels in 
the analysis ("silver channels") [|o[|3|. 

The discussion shows the strength of the analytic 
approximations, which allow to understand the 
complicated parameter interdependence. It also 
helps to optimally plan experimental setups and 
to find strategies to resolve the degeneracies. 

5. Event Rates 

The experimentally observed events must be 
compared with the theoretical expressions, which 
depend only indirectly on the above oscillation 
probabilities Every event can be classified 
by the information on the flavor of the detected 
neutrino and the type of interaction. The parti- 
cles detected in an experiment are produced by 
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neutral current (NC), inelastic charged current 
(CC) or quasi-elastic charged current (QE) inter- 
actions. The contribution to each mode depends 
on a number of factors, like detector type, the 
neutrino energy and flavour. To determine realis- 
tic event rates we compute first for each neutrino 
flavor and energy bin the number of events for 
each type of interaction in the fiducial mass of 
an ideal detector. Next the deficiencies of a real 
detector are included, like limited event recon- 
struction capabilities. The combined description 
leads to the differential event rate spectrum for 
each flavor and interaction mode as it would be 
seen by a detector which is able to separate all 
these channels. Finally different channels must 
be combined, if they can not be observed sep- 
arately. This can be due to physics, e.g., due 
to the flavor-blindness of NC interactions, or it 
can be a consequence of detector properties, e.g., 
due to charge misidentification. In order to in- 
clude backgrounds, the channels are grouped in 
an experiment specific way into pairs of signal 
and background. The considered backgrounds are 
NC-events which are misidentified as CC-events 
and CC-events identified with the wrong flavor or 
charge. For superbeams the background of CC- 
events coming from an intrinsic contamination of 
the beam is included. Finally all available signal 
channels are combined in a global analysis in or- 
der to optimally extract the physics parameters. 
The relevant channels are for a neutrino factory 
for each polarity of the beam the :^^-CC chan- 
nel (disappearance) and 77^-CC channel (appear- 
ance) event rate spectra. The backgrounds for 
these signals are NC events for all flavours and 
misidentifled v^-CC events. For superbeam ex- 
periments the signal is for each polarity of the 
beam given by the j/^-QE channel (disappear- 
ance) and z/g-CC channel (appearance). The 
backgrounds are here NC events for all flavors, 
misidentifled Vf^-CC events, and, for the Ve-CC 
channel, the I'e^CC beam contamination. 

6. The Considered LBL Setups 

The discussed sources and detectors allow dif- 
ferent LBL experiments and it is interesting to 
compare their physics potential on an equal and 



as realistic as possible footing. Studies at the 
level of probabilities are not sufficient and the 
true potential must be evaluated at the level of 
event rates as described in section ^, with realis- 
tic assumptions about the beams, detectors and 
backgrounds. We present now the results of such 
an analysis which is essentially based on refer- 
ence , where we calculate the oscillation prob- 
abilities with the exact three neutrino oscillation 
formulae in matter numerically, i.e. we use the 
approximations for the probabilities in eqs. (|^) 
and (H) only for a qualitative understanding. All 
results shown below are therefore not affected by 
approximations which were made in the deriva- 
tion of the approximate analytic oscillation for- 
mulae eqs. (|l|) and (||). Sensitivities etc. are de- 
fined by the ability to re-extract the physics pa- 
rameters from a simulation of event rates. There- 
fore event rate distributions are generated for any 
possible parameter set. Subsequently a combined 
fit to these event rate distributions is performed 
simultaneously for the appearance and disappear- 
ance channels for both polarities. This proce- 
dure uses all the available information in an opti- 
mal way. It includes spectral distributions when 
present, and it reduces to a fit of total rates when 
the total event rates are small. Adequate sta- 
tistical methods as described in are used in 
order to deal with event distributions which have 
in some regions small event rates per bin. Sys- 
tematical uncertainties are parametrized and ex- 
ternal input from geophysics is used in form of 
the detailed matter profile and its errors, which 
are included in the analysis The ability to 
re-extract in a simulation of the full experiment 
the input parameters which were used to generate 
event rate distributions is used define sensitivities 
and precision. 

It is important to include in the analysis exter- 
nal information. The discussed LBL experiments 
could in principle measure the solar Am^i and 
mixing angle 9 12. However, the precision which 
can be obtained can not compete with the ex- 
pected measurement of KamLand [Q. We in- 
clude therefore as external input the assumption 
that KamLand measures the solar parameters in 
the center of the LMA region with typical errors. 
Otherwise all unknown parameters (like the CP 
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phase) will not be constrained and are therefore 
left free, with all parameter degeneracies and er- 
ror correlations taken into account. All nuisance 
parameters are integrated out and a projection 
on the parameter of interest is performed. Alto- 
gether we are dealing with six free parameters. 

We include the beam characteristics of the 
three considered sources as well as uncertainties 
of these beam parameters, i.e. for the conven- 
tional JHF and NuMI off-axis beams uncertain- 
ties in the -background and for all beams flux 
uncertainties As detectors we con- 

sider water Cherenkov detectors, low-Z calorime- 
ters and magnetized iron detectors. For magne- 
tized iron calorimeters it is important to include 
realistic threshold effects. We use a linear rise 
of the efficiency between 4 GeV and 20 GeV and 
we study the sensitivity to the threshold position. 
We do not include liquid Argon TPCs in our anal- 
ysis, but they would certainly be an important 
detector if this technology will work. The con- 
sidered beams and detectors allow now different 
interesting combinations which are listed in ta- 
ble 1^. JHF-SK is the planned combination of the 
existing SuperKamiokande detector and the JHF 
beam, while JHF-HK is the combination of an 
upgraded JHF beam with the proposed Hyper- 
Kamiokande detector. With typical parameters, 
JHF-HK is altogether about 95 times more inte- 
grated luminosity than JHF-SK, and we assume 
that it operates partly with the anti-neutrino 
beam. Water Cherenkov detectors are ideal for 
the JHF beam, since charged current quasi elas- 
tic scattering is dominating. A low-Z calorimeter 
is proposed for the NuMI off-axis beam, which is 
better here, since the energy is higher and there 
is already a considerable contribution of inelas- 
tic charged current interactions. NuFact-I is an 
initial neutrino factory, while NuFact-II is a fully 
developed machine, with 42 times the luminosity 
of NuFact-I 0,^9],^ . Deep inelastic scattering 
dominates for these even higher energies and mag- 
netized iron detectors are therefore considered in 
combination with neutrino factories. 



acronym 


detector 


L 


L/ -Bpcak 


JHF-SK 


water Cherenkov 


295 


378 


NuMI 


low-Z 


735 


337 


NuFact-I 


10 kt mag. iron 


3000 


90 


JHF-HK 


water Cherenkov 


735 


295 


NuFact-II 


40 kt mag. iron 


3000 


90 



Table 1 

The considered combinations of beams and detec- 
tors and their acronyms. 



7. Results 

A realistic analysis of future LBL experiments 
requires a number of different aspects to be taken 
into account. It should be clear from the discus- 
sion above that it is not sufficient to quote lim- 
its which are based on oscillation probabilities or 
merely on the statistics of a single channel with- 
out backgrounds or systematics. Depending on 
the position in the space of physics parameters 
the degeneracies or correlations, the backgrounds, 
the systematics or statistics may be the limiting 
factor. A reliable comparative study of the dis- 
cussed LBL setups requires therefore a detailed 
analysis of the six dimensional parameter space, 
which includes all these effects on the same foot- 
ing @. ^ 

There is excellent precision for the leading os- 
cillation parameters Am|j and sin'^2023, which 
will not be further discussed here. The more in- 
teresting sensitivity to the sub-leading parameter 
sin^ 2013 depends on what will be found for Am^i 



and Amjj. Assuming that the leading para- 
meters are measured to be Ato§i = 3 • 10"'^ eV^, 
sin^ 2^23 = 0.8 and that KamLand measures the 
solar parameters at the current best fit point of 
the LMA region, i.e. Am|i = 6 • 10"^ eV^ and 
sin 2^12 = 0.91, we can make a comparison of the 
sin^ 2013 sensitivity limit for the different setups. 
The result is shown in fig. ^j. The individual con- 
tributions of different sources of uncertainties are 
shown for every experiment and the left edge of 
every band in fig. ^ corresponds to the sensitivity 
limit which would be obtained purely on statisti- 
cal grounds. This limit is successively reduced by 
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Sensitivity to sin 2^13 



I Systematic 
Correlation 
Degeneracy 



10-' 



10" 



JHF-SK 

NuMI 

JHF-HK 

NuFact-I 

NuFact-II 



10"'' 10"' 
sin^2ei3 



10" 



10" 



Sensitivity to the sign of Ami 



I Systematic 
I Cotitlation 



10-' 



lo- 



10-* IQ-' 



JHF-SK 
^Ni.)MI 
JHF-HK 
NuFact-I 
NuFact-II 
10"^ 10"' 



Figure 3. The sin^ 26'i3 sensitivity for all setups 
defined in section at the 90% confidence level 
for Aml^ ^ 3 • IQ-^eV^ and sin^ 26*23 = 0.8. The 
plot shows the deterioration of the sensitivity lim- 
its as the different error sources are successively 
switched on. The left edge of the bars is the sen- 
sitivity statistical limit. This limit gets reduced 
as systematical, correlational and degeneracy er- 
rors are switched on. The right edge is the final 
sensitivity limit |Q. 



adding the systematical uncertainties of each ex- 
periment, the correlational errors and finally the 
degeneracy errors. The right edge of each band 
constitutes the final error for the experiment un- 
der consideration. It is interesting to see how the 
errors of the different setups are composed. There 
are different sensitivity reductions due to system- 
atical errors, correlations and degeneracies. The 
largest sensitivity loss due to correlations and de- 
generacies occurs for NuFact-II, which is mostly 
a consequence of the uncertainty of Am^i , which 
translates into an a uncertainty, and which domi- 
nates the appearance probability eq. (|^) for small 
sin^ 2013. Note that it is in principle possible to 
combine different experiments. If done correctly, 
this allows to eliminate part or all of the correla- 
tional and degeneracy errors [|8|. 

Another challenge of future LBL experiments 
is to measure signlArn^i) via matter effects and 



Figure 4. The sin^ 26*13 sensitivity region to 
si(7rt(Am|j^) for the setups defined in section |^. 
The left edge of the bars are the statistical sen- 
sitivity limits which are successively reduced by 
systematical, correlational and degeneracy errors. 
The right edge of the bars is the final limit. 



the sensitivity which can be obtained for the se- 
tups under discussion is shown in fig. ^. Taking 
all correlational and degeneracy errors into ac- 
count we can see that it is very hard to de- 
termine sign{Am^i) with the considered super- 
beam setups. The main problem is the degen- 
eracy with S, which allows always the reversed 
sign{Am^i) for another CP phase. Note, how- 
ever, that the situation can in principle be im- 
proved if different superbeam experiments were 
combined such that this degeneracy error could 
be removed. Neutrino factories perform consider- 
ably better on sign{Am'^i), particularly for larger 
baselines. Combination strategies would again 
lead to further improvements. 

Coherent forward scattering of neutrinos and 
the corresponding MSW matter effects are so far 
experimentally untested. It is therefore very im- 
portant to realize that matter effects will not only 
be useful to extract sign{Amli), but that they al- 
low also detailed tests of coherent forward scatter- 
ing of neutrinos. This has been studied in detail 
in ||,||4|,||. 

The Holy Grail of LBL experiments is the mea- 
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Sensitivity to CP- Violation at 6cp = +ir/2 



- LlxlO"-^ JHF-HK 




10"' lO-"" 10"* 10"* 10"' 10"^ 10-1 jqO 
sin^2ei3 



Figure 5. The sin^ 26*13 sensitivity range for CP 
violation of the considered setups at 90% confi- 
dence level and for different Am^i values. The 
upper row corresponds to the lower bound of 
Amli = 1.1 X 10~^ eV^, the bottom row to 
the upper bound Am|i = 4.7 x 10""* eV^ and 
the middle row to the best LMA fit, Atoji = 
3.7 X 10~^ eV^. Cases which do not have CP sen- 
sitivity are omitted from this plot. The chosen 
parameters are S = +n/2, Am§i = 3 • 10"'^ eV^, 
sin^ 26*23 = 0.8, and a solar mixing angle cor- 
responding to the current best fit in the LMA 
regime §. 



surement of leptonic CP violation. The sin^ 2^13 
sensitivity range for measurable CP violation is 
shown in fig. ^ for ^ = 7r/2 for the different se- 
tups and for different values of Arn^i- It can be 
seen that measurements of CP violation are in 
principle feasible both with high luminosity su- 
perbeams as well as advanced neutrino factories. 
However, the sensitivity depends in a crucial way 
on Ato|i. For a low value Am^i = 1.1 10"^ eV^, 
the sensitivity is almost completely lost, while the 
situation would be very promising for the largest 
considered value Arriji = 4.7 10""* eV^. For a 
measurement of leptonic CP violation it would 
therefore be extremely exciting and promising if 
KamLand would find Am^i on the high side of 
the LMA solution (the so-called HLMA case). 
The sensitivities shown in fig. depend on the 



choice for S. The value which was used here was 
6 — 11/2 and and the sensitivities become become 
worse for small CP phases close to zero or tt. 

8. Conclusions 

Future long baseline neutrino oscillation exper- 
iments will lead to precision neutrino physics. A 
basic fact which makes this possible is that the 
atmospheric mass splitting Arri^i — Arn^^^ leads 
for typical neutrino energies iJ^ ~ 1 — 100 GeV 
to oscillation baselines in the range 100 km to 
10000 km. Beams have the advantage, that unlike 
the sun or the atmosphere, they can be controlled 
very precisely. Combined with equally precise 
detectors and an adequate oscillation framework 
(including three neutrinos and matter effects) 
must be used. There exist different interesting 
sources for long baseline oscillation experiments, 
like reactors or /3-beams, but we restricted the 
discussion here to superbeams and neutrino fac- 
tories. We presented the issues which enter into 
realistic assessments of the potential of such ex- 
periments. The discussed experiments turned out 
to be very promising and they lead to very pre- 
cise measurements of the leading oscillation para- 
meters Am§i and sin^2^23- We discussed in de- 
tail how the different setups lead to very interest- 
ing measurements or limits 613 and 6. It will also 
be possible to perform impressive tests of Earth 
matter effects, allowing to extract sign{Arri^i) . 
The discussed setups have an increasing poten- 
tial and increasing technological challenges, but 
it seems possible to built them in stages. The 
shown results are valid for each individual setup 
and future results should of course be included in 
the analysis. This would be especially important 
if more LBL experiments were built and depend- 
ing on previous results there exist different opti- 
mization strategies. In the short run the expected 
results from KamLand are extremely important 
and will have considerable impact. First it will 
become clear if Am|i lies in the LMA regime, 
which is very important since realistically, CP vi- 
olation can only be measured then. Within the 
LMA solution it is also very important if Atoj^ 
lies close to the current best fit, on the high or on 
the low side. A value of Am2i on the high side 
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(HLMA) would be ideal, since it would guaran- 
tee an extremely promising LBL program with a 
chance to see leptonic CP violation already with 
the JHF beam in the next decade. 
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